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Abstract 



The action and topological charge are used to determine the relative rates of 
standard cooling and smearing algorithms in pure SU(3)-color gauge theory. 
We consider representative gauge field configurations on 16'^ x 32 lattices at 
(3 = 5.70 and 24^ x 36 lattices at /? = 6.00. We find the relative rate of 
variation in the action and topological charge under various algorithms may 
be succinctly described in terms of simple formulae. The results are in accord 
with recent suggestions from fat-link perturbation theory. 
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I. INTRODUCTION 



Cooling is a well established method for locally suppressing quantum fluctuations 
in gauge field configurations. The locality of the method allows topologically nontrivial 
field configurations to survive numerous iterations of the cooling algorithm. Application of 
such algorithms are central to removing short-distance fluctuations responsible for causing 
renormalization constants to significantly deviate from one as is the case for the topological 
charge operator 0]. 

An alternative yet somewhat related approach is the application of APE smearing 
to the gauge field configurations. Because the algorithm is applied simultaneously to all 
link variables, without any annealing in the process of a lattice update, the application of 
the algorithm may be viewed as an introduction of higher dimension irrelevant operators 
designed to improve the operator based on unsmeared links. This approach also leads to the 
desirable effect of renormalization constants approaching one. The net effect may be viewed 
as introducing a form factor to gluon vertices 0] that suppresses large interactions and, 
as a result, suppresses lattice artifacts at the scale of the cutoff. 

The focus of this paper is to calibrate the relative rates of cooling and smearing as mea- 
sured by the action and the topological charge. In doing so we aim to establish relationships 
between various algorithms and their internal parameters. It is well known that errors in 
the standard Wilson action eventually destroy (anti)instanton configurations and this effect 
is represented in the topological charge as sharp transitions from one integer to another. We 
find the timing analysis of these topological charge defects to be central to understanding 
the behavior of the various algorithms. 

To bridge the gap between cooling and smearing and discover the origin of variations 
between these algorithms, we present a new algorithm for smearing the gauge fields which 
we refer to as Annealed U Smearing or AUS smearing. In addition to comparing cooling and 
smearing, we also calibrate the effects of different smearing fractions, describing the balance 
between the unsmeared link and the smeared link. 

There have been numerous studies of the dependence of vacuum matrix elements on the 
number of smoothing sweeps. Of particular note is the observation of cooling as a diffusive 
procedure 0, such that the distance d affected by cooling is related to the square-root of 
the number of sweeps over the lattice, n. This random walk ansatz was further investigated 
by introducing the lattice spacing a, such that d = a^pn takes on a physical length p. With 
this simple formula, it is possible to effectively relate results from lattices of different (5 (a) 
and different smoothing sweeps n. Extensive work on removing the dependence of particular 
quantities on the number of smoothing sweeps was carried out in Ref. 0. There, highly 
improved actions are fine tuned to stabilize the instanton action and size under cooling. 
Moreover, a modification of the updating algorithm [Q allows one to define a physical length 
scale beyond which the physics is largely unaltered. Closer to the focus of this investigation, 
some exploratory comments favoring the possibility of calibrating under-relaxed cooling were 
made in Ref. 0. We will lend further weight to this possibility in the following analysis. 

This paper is divided as follows: Section briefiy describes the gauge action and sum- 
marizes the numerical simulations examined in this study. Section |ni| presents the relaxation 
algorithms and topological charge definition on the lattice. Sections [IV A| and [IV B| present 
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the numerical results for the action and topological charge trajectories respectively. A sum- 
mary of our findings is provided in Section 0. 



II. LATTICE GAUGE ACTION FOR SUc{3) 
A. Gauge Action. 

The Wilson action is defined as, 

Sg = PY.J2 \neTT{l - U,Ax)), (1) 

where the operator U^^{x) is the standard plaquette, 

U^^uix) = U,{x) U,{x + fi) Ulix + v) Ulix) . (2) 

Gauge configurations are generated using the Cabbibo-Marinari pseudo-heat-bath 
algorithm with three diagonal SUc{2) subgroups. Simulations are performed using a parallel 
algorithm on a Thinking Machines Corporations (TMC) CM-5 with appropriate link parti- 
tioning. For the Wilson action we partition the link variables in the standard checkerboard 
fashion. 

Configurations are generated on a 16^ x 32 lattice at /5 = 5.70 and a 24^ x 36 lattice at 
(3 = 6.00. Configurations are selected after 5000 thermalization sweeps from a cold start, 
and every 1000 sweeps thereafter. Lattice parameters are summarized in Table |l[ When 
required, we use the plaquette definition of the mean link variable. 



III. COOLING AND SMEARING ALGORITHMS 

Relaxation techniques have been the subject of a large amount of study in lattice QCD 
over the past few years [|lT|. Numerous algorithms have been employed to suppress the 



high-frequency components of the gauge field configurations in investigating their semi- 
classical content. In the following three subsections we consider three different algorithms 
that perform this task. 

Perhaps it should be noted that cooling algorithms have been designed to render in- 
stantons stable over many hundreds of sweeps. For example, De Forcrand et al. have 
established a five-loop improved gauge action that achieves this. However, our interest 
here is in using the defects in the topological charge evolution under standard cooling and 
smearing algorithms based on the Wilson gauge action to calibrate these algorithms. 

TABLE L Parameters of the numerical simulations. 



Action 


Volume 


^Thcrm 


^Samp 


/? 


a (fm) 




Physical Volume (fm) 


Wilson 


16^ X 32 


5000 


1000 


5.70 


0.18 


0.86085 


2.88^ X 5.76 


Wilson 


X 36 


5000 


1000 


6.00 


0.10 


0.87777 


2.43 X 3.6 
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A. Cooling 



Standard cooling minimizes the action locally at each link update. The preferred algo- 
rithm is based on the Cabbibo-Marinari |jTO| pseudo-heat-bath algorithm for constructing 



SUc{S)-co\or gauge configurations. A brief summary of a link update follows. 

At the SUc{2) level the algorithm is transparent. An element of SUc{2) may be param- 
eterized as, U = Oq/ + ia ■ a, where a is real and = 1. Let be one of the six staples 
associated with creating the plaquette associated with a link U^. 

tj, = U,{x + fi) Ulix + v) Ul{x) . (3) 

We define 



Y,U^ = kU where U G SUc{2) and fc^ = det ^ f/„ . (4) 




The feature of the sum of SUc{2) elements being proportional to an SUc{2) element is central 
to the algorithm. The local SUc{2) action is proportional to 

7^eTr(l-f/F), (5) 

and is locally minimized when TZe Tt{UU) is maximized, i.e. when 

7^eTr([/^7) =7^eTr(/). (6) 

which requires the link to be updated as 



U-^U' = =W= ^ ^ . (7) 

rC 

At the SUc{3) level, one successively applies this algorithm to various SUc{2) subgroups of 
SUc{3), with SUc{2) subgroups selected to cover the SUc{3) gauge group. 

We explored cooling gauge field configurations using two diagonal SUc{2) subgroups. 
The cooling rate is slow and the action density is not smooth, even after 50 sweeps over 
the lattice. Addition of the third diagonal SUc{2) subgroup provides acceptably fast cooling 
and a smooth action density. These SUc{2) subgroups are acting as a covering group for the 
SUc{3) gauge group. One would therefore expect that repeated updates of SUc(2) subgroups 
will better approximate the ultimate SUc{'i) link. 



B. Smearing 

Here we consider two algorithms for smearing the gauge links. APE smearing is now 
well established. Our alternate algorithm AUS smearing is designed to remove instabilities of 
the APE algorithm and provide an intermediate algorithm sharing features of both smearing 
and cooling. 
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1. APE Smearing 



APE smearing is a gauge equivariant prescription for averaging a link UJx) with 



its nearest neighbors f/^(a; + z>), p ^ fi. The hnear combination takes the form: 



X 



o 



(8) 



where S(x;/i) = [Y,l, 



X 



is the sum of the six staples defined in Eq. 



The 



parameter a represents the smearing fraction. The algorithm is constructed as follows: (i) 
calculate the staples, S(a;; /i); (ii) calculate the new link variable U'^{x) given by Eq. (|) and 
then reunitarize U'^{x); (iii) once we have performed these steps for every link, the smeared 
U'^{x) are mapped into the original U^{x). This defines a single APE smearing sweep which 
can then be repeated. 

The celebrated feature of APE smearing is that it can be realized as higher-dimension 
operators that might appear in a fermion action for example. Indeed, "fat link" actions 
based on the APE algorithm are excellent candidates for an efficient action with improved 
chiral properties [|13|. The parameter space of fat link actions is described by the number of 
smearing sweeps riape and the smearing coefficient a. It is the purpose of this investigation 
to explore the possible reduction of the dimension of this parameter space from two to one. 



2. A US Smearing 

The Annealed U Smearing (AUS) algorithm is similar to APE smearing in that we take a 
linear combination of the original link and the associated staples, /i). However in AUS 
smearing we take what was a single APE smearing sweep over all the links on the lattice 
and divide it up into four partial sweeps based on the direction of the links. One partial 
sweep corresponds to an update of all links oriented in one of the four Cartesian directions 
denoted /z = 1,..,4. Hence in a partial sweep we calculate the reuniterized U'^{x) for all 
/i-oriented links and update all of these links [U'^{x) — > Ufj,{x)] at the end of each partial 
sweep. 

Thus the difference between AUS smearing and APE smearing is that in APE smearing 
no links are updated until all four partial sweeps are completed, while in AUS smearing, 
updated smeared information is cycled into the calculation of the next link direction, a 
process commonly referred to as annealing. In a sense AUS smearing is between cooling 
and APE smearing in that cooling updates one link at a time, AUS smearing updates one 
Cartesian direction at a time, and APE smearing updates the whole lattice at the same 
time. 

As for APE smearing, the reunitarization of the SUc{3) matrix is done using the standard 
row by row orthonormalization procedure: begin by normalizing the first row; then update 
the second row by row2 = row2 — {row2 ■ rowl)rowl; normalize the second row; finally set 
row3 equal to the cross product of rowl and row2. 

For SUc{2) gauge theory, cooling and smearing are identical for the case of a = 1 when 
the links are updated one at a time. The AUS algorithm changes the level of annealing 
and provides the opportunity to alter the degree of smoothing. At the same time the 
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AUS algorithm preserves the gauge equivariance of the APE algorithm for SUc{S) gauge 
theory 
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As we shall see, it also removes the instability of the APE algorithm at large 
smearing fractions a. This latter feature may provide a more efficient method for finding the 
fundamental modular region of Landau Gauge in trying to understand the effects of Gribov 



copies in gauge dependent quantities such as the gluon-propagator |]T^ 



C. Topological Charge 

We construct the lattice topological charge density operator analogous to the standard 
Wilson action via the clover definition of F^^^. The topological charge Ql and the topological 
charge density qh{x) are defined as follows 



327r- 



(9) 



where the field strength tensor is 



a^gF. 



3Tr(0^.(x) 



(10) 



and 0^y{x) is 



0^,{x] 



= U^{x) Uy{x + /t) Ul{x + z>) Ul{x) 

+ Uy{x) Ul{x + v - fi) Ul{x - fi) Uf,{x - fi) 

+ U^{x — fi) Ul{x — fi — 0) U^{x — fi — 0) Uy{x — z>) 

+ Ul{x - u) U^{x - 0) U^{x + fi-0) UUx) . 



(11) 



Lattice operators possess a multiplicative lattice renormahzation factor, Ql = Zq{(3)Q 
which relates the lattice quantity to the continuum quantity Q. Perturbative calculations 
indicate Zq{(3) ^ 1 — 5.451//3 + 0{1/ j3'^) 0]. This large renormahzation causes a problem 
when one is working at /5 ~ 6.0, as Zq{(3) <ti 1. This implies that the topological charge 
is almost impossible to calculate directly. However when one applies cooling or smearing 
techniques to remove the problem of the short range quantum fiuctuations giving rise to 
Zq{P) <^ 1, one can resolve near integer topological charge. One can apply the operator 
Ql to cooled configurations or smeared configurations. The latter case may be regarded as 
employing an improved operator in which the smeared links are understood to give rise to 
additional higher dimension irrelevant operators designed to provide a smooth approach to 
the continuum limit. 



IV. NUMERICAL SIMULATIONS 

We analyzed two sets of gauge field configurations generated using the Cabbibo- 
Marinari [|1^] pseudo-heat-bath algorithm with three diagonal SUc{2) subgroups as described 
in Table |. Analysis of a few configurations proves to be sufficient to resolve the nature of 
the algorithms in question. The two sets are composed of ten 16^ x 32 configurations and 
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five 24^ X 36 configurations. For each configuration we separately performed 200 sweeps of 
cooling, 200 sweeps of APE smearing at four values of the smearing fraction and 200 sweeps 
of AUS smearing at six values of the smearing fraction. 

For clarity, we define the number of times an algorithm is applied to the entire lattice 
as ric, ?T,ape(a) and naus(tt) for cooling, APE smearing and AUS smearing respectively. We 
monitor both the total action normalized to the single instanton action 5*0 = Svr^/gf^ and 
the topological charge Qh{x) of Eq. and observe their evolution as a function of the 
appropriate sweep variable and smearing fraction a. 

For APE smearing we consider four different values for the smearing fraction a includ- 
ing 0.30, 0.45, 0.55, and 0.70. Larger smearing fractions reveal an instability in the APE 
algorithm where the links are rendered to noise. 

The origin of this instability is easily understood in fat-link perturbation theory |^ where 
the smeared vector potential after n APE smearing steps is given by 



reflecting the transverse nature of APE smearing. Here 



4"^?), (12) 



g. = -sinf^) (13) 



and 

j\q) = l- f. (14) 



For /(g) to act as a form factor at each vertex in perturbation theory over the entire Brillouin 
zone 

TT TT 

--<q,<-. (15) 
a a 

one requires —1 < f{q) < 1 which constrains a to the range < a < 3/4. 

The annealing process in AUS smearing removes this instability. Hence, the parameter 
set for AUS smearing consists of the APE set plus an extra two, a = 0.85 and 1.00. 



A. The Action Analysis. 

1. APE and AUS Smearing Calibration 

The action normalized to the single instanton action S/ Sq can provide some insight into 
the number of instantons left in the lattice as a function of the sweep variable for each 
algorithm. However, our main concern is the relative rate at which the algorithms perform. 
In Fig. H, we show S/Sq as a function of cooling sweep on the 24^ x 36 lattice for five 
different configurations. The close proximity of the five curves is typical of the configuration 
dependence of the normalized action S/Sq. 
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FIG. 1. The ratio S/Sq as a function of cooling sweeps ric for five configurations on the 24^ x 36 
lattice at P = 6.00. The single instanton action is 5*0 = Sir^/g'^. 
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FIG. 2. The ratio S/Sq as a function of APE smearing sweeps nape (a) for a configuration on a 
24^ X 36 lattice at /? = 6.00. Each curve has an associated smearing fraction a. 
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lO' F 




50 100 150 200 

n (a) 

aus ^ ^ 

FIG. 3. The ratio S/Sq as a function of AUS smearing sweep n.^^s{o.) for a configuration on a 
24^ X 36 lattice at P = 6.00. Each curve has an associated smearing fraction a. 

Figures |^ and |^ report results for APE and AUS smearing respectively. Here we fo- 
cus on one of the five configurations, noting that similar results are found for the other 
configurations. Each curve corresponds to a different value of the smearing fraction a. 

A similar analysis of the 16^ x 32 lattice at /3 = 5.70 is also performed yielding analo- 
gous results. In fact, taking the physical volumes of the two lattices into account reveals 
qualitatively similar action densities after 200 sweeps. 

Note that in Fig. |^, the curve associated with the smearing parameter a = 1.00, crosses 
over the one generated at a = 0.85, when the sweep number, riausl^), is approximately 
40 sweeps. Thus a = 1 is not the most efficient smearing fraction for untouched gauge 
configurations. However, as the configurations become smooth, a = 1.00 becomes the most 
efficient. 

To calibrate the rate at which the algorithms reduce the action, we record the number of 
sweeps n{a) required for the smeared action to cross various thresholds St- This is repeated 
for each of the smearing fractions a under consideration. In establishing the relative a 
dependence for the number of sweeps n{a), we first consider a simple linear relation between 
the number of sweeps required to cross 5*^ at one a compared to another a', i.e. 

n{a') = Co + ci n{a) . (16) 

Anticipating that Cq will be small if not zero, we divide both sides of this equation by n{a) 
and plot as a function of n{a). Deviations from a horizontal line will indicate failings of our 
linear assumption. 

Fig. ^ displays results for a' fixed to 0.55 for the APE smearing algorithm and Fig. |^ 
displays analogous results for AUS smearing. We omit thresholds that result in n(0.55) < 
10 as these points will have integer discretization errors exceeding 10%. For a = 0.70 
discretization errors the order of 10% are clearly visible in both figures. 
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FIG. 4. The ratio nape(0.55)/nape(Q;) versus nape(a) for numerous threshold actions on the 
24^ X 36 lattice at /3 = 6.00. Prom top down the data points correspond to a = 0.70, 0.55, 0.45 
and 0.30. 
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FIG. 5. The ratio naus(0.55)/naus(Q;) versus naus(«) for numerous threshold actions on the 
24^ X 36 lattice at /? = 6.00. Prom top down the data points correspond to a = 1.00, 0.85, 0.70, 
0.55, 0.45, 0.30. 
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a 

FIG. 6. Illustration of the a dependence of ci = (n(a' = 0.55)/n(a)) for APE smearing. The 
solid line fit to the data indicates ci = 1.838 a — 0.011 whereas the dashed line, constrained to pass 
through the origin, provides a slope of 1.818. 



For the smearing fraction a < 0.85 both plots show little dependence of n{a' = 
0.55) /n{a) on n{a). This supports our simple ansatz of Eq. (p!6| ) and indicates cq is in- 
deed small as one would expect. 

The non-linear behavior for a = 1.00 in AUS smearing, in Fig. |^, reflects the cross 
over in Fig. ^ for a = 1.00. For unsmeared configurations, a = 1.00 appears to be too 
large. Further study in SUc{2) is required to resolve whether the origin of the smearing 
inefficiency lies in the sum of staples lying too far outside the SUd^S) gauge group for useful 
reunitarization, or whether the annealing of the links in AUS smearing is insufficient relative 
to cooling. 

To determine the a dependence of ci, we plot ci = {n{a' = 0.55)/n(a)) as a function of 
a. Here the angular brackets denote averaging over all the threshold values considered; i.e. 
averaging over data in the horizontal lines of Figs. ^ and ^. 

Fig. 1^ for APE smearing and Fig. |^ for AUS smearing indicate the relationship between 
Ci and a is linear with zero intercept to an excellent approximation. Indeed, ignoring the 
point at a = 1 for AUS smearing, one finds the same coefficients for the a dependence of 
APE and AUS smearing. When the fits are constrained to pass through the origin, one finds 
a slope of 1.818 which is the inverse of a' = 0.55. Hence we reach the conclusion that 



'■ape V , ^ 



7 (17) 



This analysis based on the action suggests that a preferred value for a does not really 
exist. In fact it has been recently suggested that one should anticipate some latitude in the 
values for n{a) and a that give rise to effective fat-link actions [Q. What we have done here 
is established a relationship between n{a) and a, thus reducing what was potentially a two 
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a 

FIG. 7. Illustration of the a dependence of ci = (n(a' = 0.55)/n(a)) for AUS smearing. Fits 
to the data exclude the point at a = 1. The solid line fit to the data indicates ci = 1.857 a — 0.021 
whereas the dashed line, constrained to pass through the origin, provides a slope of 1.818. 



dimensional parameter space to a one dimensional space. This conclusion will be further 
supported by the topological charge analysis below. 
To summarize these finding we plot the ratios 



a' nape (a') , a' n^^^ia' 



and — — , (18) 

designed to equal 1 in Figs. || and |^ for APE and AUS smearing respectively. Figs. ^ and 



11 report the final results of a similar analysis for APE and AUS smearing respectively at 



P = 5.7. Here the a = 1.0 results are omitted from the AUS smearing results for clarity. 



2. Cooling Calibration 

Here we repeat the previous analysis, this time comparing cooling with APE and AUS 
smearing. We make the same linear ansatz for the relationship and plot ric/na 



^ape(tt) for APE smearing in Fig. |T2|. Fig. |T3| reports the ratio nc/riaus(«) versus naus(«) for 



AUS smearing. At small numbers of smearing sweeps, large integer discretization errors the 
order of 25% are present, as Uc is as small as 4. With this in mind, we see an independence 
of the ratio on the amount of cooling/smearing over a wide range of smearing sweeps for 
both APE and AUS smearing. This supports a linear relation between the two algorithms. 
Averaging the results provides 

^'^ 0.330 and - = 0.340 , (19) 
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(0.55) 

or more generally 
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FIG. 8. Illustration of the degree to which the relations of ( [T7| ) are satisfied for the action under 
APE smearing. Here a < a' and a' = 0.30, 0.45, 0.55, and 0.70. Data are from the 24^ x 36 lattice 
at /? = 6.0. 
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FIG. 9. Illustration of the degree to which the relations of (17) are satisfied for the action under 
AUS smearing. Here a < a' and a' = 0.30, 0.45, 0.55, 0.70 and 0.85. Data are from the 24^ x 36 
lattice at (3 = 6.0. 



13 



2.0 



1.8 - 

1.6 - 

1.4 - 

1.2 - 

1.0 - 

0.8 - 



gf' 0.6 - 

^ 0.4 - 

0.2 - 

0.0 I ^ ^ ^ ^ ^ ^ ^ ^ ^ — 

20 40 60 80 100 120 140 160 180 200 

n (a) 

FIG. 10. Illustration of the degree to which the relations of ( p^ are satisfied for the action 
under APE smearing. Here a < a' and a' = 0.30, 0.45, 0.55, and 0.70. Data are from the 16^ x 32 
lattice at /3 = 5.7. 
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FIG. 11. Illustration of the degree to which the relations of ( p^ are satisfied for the action 
under AUS smearing. Here a < a' and a' = 0.30, 0.45, 0.55, 0.70 and 0.85. Data are from the 
16^ X 32 lattice at /? = 5.7. 
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FIG. 12. Ratio of cooling to APE smearing sweeps as a function of the number of APE smearing 
sweeps. From top down a = 0.70 0.55, 0.45 and 0.30. 
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FIG. 13. Ratio of cooling to AUS smearing sweeps as a function of the number of AUS smearing 
sweeps. From top down a = 1.00, 0.85 0.70 0.55, 0.45 and 0.30. 
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nc ~ 0.600 a nape («) and nc — 0.618 anaus(«) • (20) 

Hence we see that cooling is much more efficient at smoothing than the smearing algorithms 
requiring roughly half the number of sweeps for a given product of n and a. 

Equating the equations of Eq. ( pOD provides the following relation between APE and 
AUS smearing 



1.03«'naus(«'), (21) 



summarizing the near equivalence of the two smearing algorithms. It is important to recall 
that the annealing of the AUS smearing algorithm removes the instability of APE smearing 
encountered at large smearing fractions, a. Thus AUS smearing offers a more stable gauge 
equivariant smoothing of gauge fields. It offers faster smoothing due to its ability to handle 
larger smearing fractions. This algorithm may be of use in studying Gribov ambiguities in 
Landau gauge fixing |]12 . 



A similar analysis of the 16^ x 32 lattice at /? = 5.7 provides 

nc ~ 0.572 a riapc (a) and nc ~ 0.604anaus(a) , (22) 

with 

a^ape(a) - 1.06a'naus(a') • (23) 

Here the change in the coefficient relating APE smearing to cooling appears to be propor- 
tional to (3. 



B. The Topological Charge Density Analysis. 

Typical evolution curves for the lattice topological charge operator of Eq. (0) are shown 



in Figs. |14|, |T5| and |T6| for APE smearing, AUS smearing and cooling respectively. These 
data are obtained from a typical gauge configuration on our 24'^ x 36 lattice at /3 = 6.00. 
The configuration used here is the same representative configuration illustrated in Figs. |^ 
and ^ of the action analysis. 

The main feature of these figures is that the smearing/cooling algorithms produce a 
similar trajectory for the topological charge. For most cases only the rate at which the 
trajectory evolves changes. In these cases, one can use these trajectories as another way in 
which to calibrate the rates of the algorithms. 

After a few sweeps, the topological charge converges to a near integer value with an error 
of about 10%, typical of clover definitions of Q at /3 = 6.0. The standard Wilson action is 
known to lose (anti)instantons during cooling due to 0{a?) errors in the action which act to 
tunnel through the single instanton action bound. Given the intimate relationship between 



cooling and smearing discussed in Section III it is not surprising to see similar behavior in 



the topological charge trajectories of the smearing algorithms considered here. The sharp 
transitions from one integer to another indicate the loss of an (anti)instanton. 

These curves look quite different for other configurations. Fig. ^ displays trajectories 
for cooling on five different gauge configurations. However, the feature of similar trajec- 
tories for the various cooling/smearing algorithms remains at /3 = 6.0. To calibrate the 
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FIG. 14. Typical evolution curve of the lattice topological charge operator as a function of the 
number of sweeps for APE smearing. Data are from the 24^ x 36 lattice at /3 = 6.00. Each curve 
corresponds to a particular a value and indicated. 
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FIG. 15. Typical evolution curve of the lattice topological charge operator under AUS smearing. 
Data are from the 24^ x 36 lattice at /? = 6.00. Each curve corresponds to a particular a value and 
indicated. 
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FIG. 16. Typical evolution curve of the lattice topological charge operator as a function of the 
number of sweeps for cooling with three diagonal SU(2) subgroups on a 24^ X 36 lattice at /3 = 6.00. 



cooling/smearing algorithms, we select thresholds at topological charge values where the 
trajectories are making sharp transitions from one near integer to another. Table |I| summa- 
rizes the thresholds selected and the number of sweeps required to pass though the various 
thresholds. Note that the second configuration, C2, data entry in Table fl| corresponds to 
the sampling of the curves illustrated in Figs. 14, 15 and 16. 

At /5 = 5.7 we did not find analogous trajectories, suggesting that the coarse lattice 
spacing and larger errors in the action prevent one from reproducing a similar smoothed 
gauge configuration using different algorithms. In this case the two-dimensional aspect of 
the smearing parameter space remains for studies of the topological sector. That this might 
be the case is hinted at in Fig. 11 in which the ratio of smearing results for /? = 5.7 lattices 
is not as closely constrained to one as for the (3 = 6.0 results in Figs. |^. 

As in the action analysis, we report the ratio of a = 0.55 results to other a value results 
within APE and AUS smearing. Figs. ^ and |1^ summarize the results for APE and AUS 
smearing respectively. Again we see an enhanced spread of points at small numbers of 
smearing sweeps due to integer discretization errors. 



Data for large numbers of sweeps at small and large a values are absent in Fig. 19. The 



absence of points for small a values is simply due to the thresholds not being crossed within 
200 smearing sweeps. The absence of points for large a values reflects the divergence of the 
topological charge evolution from the most common trajectory among the algorithms. This 
divergence is also apparent in Fig. ^ where the a = 1.0 AUS smearing results fail to satisfy 
the linear ansatz. 

While the data from the topological charge evolution is much more sparse, one can 
see reasonable horizontal bands forming supporting a dominant linear relationship between 
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FIG. 17. The trajectories of the lattice topological charge density operator as a function of the 
number of sweeps for cooling on typical 24^ x 36, (3 = 6.00 configurations. 
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FIG. 18. The ratio nape(0.55)/nape(a) versus nape(a) for from the data of Table || extracted 
from the 24'^ x 36 configurations at /? = 6.00. From the top down, the horizontal sets of points 
correspond to a = 0.70, 0.55, 0.45 and 0.30. 
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TABLE II. Summary of the number of sweeps required to pass through various topological 

charge thresholds. The selection of the thresholds is described in the text. Smearing fractions a 
arc indicated in the table headings. Omissions in the table indicate either the threshold was not 
met within 200 sweeps or that the trajectory diverged from the most common trajectory among 
the algorithms. 
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FIG. 19. The ratio naus(0.55)/naus(o) versus naus(a) for Ql from the data of Table ||extracted 
from the 24^ x 36 configurations at /3 = 6.00. From the bottom up, the symbols correspond to 
a = 0.30, 0.45, 0.55, 0.70, 0.85 and 1.00. 



various a values. Averages of the bands are reported in Table |T| along with previous results 
from the action analysis. With the exception of the a = 1.0 AUS smearing results, the 
agreement is remarkable, leading to the same conclusions of the action analysis summarized 
in Eq. (^). 

In calibrating cooling via the topological charge defects, we report the ratio nc/nape(0.55) 
as a function of nape(a) and ric/ ?T-aus(0.55) as a function of naus(«) for APE and AUS smearing 
respectively. The results are shown in Figs. ^ and |21]. 

The data is too poor to determine anything beyond a linear relation between ric and 
nape{(y) or naus(«)- Averaging the results provides 

-£—- = 0.30(3) and ^ = 0.35(3), (24) 

nape(0.55) naus(0.55) 

in agreement with the earlier action based results of Eq. (p!9|). 

As a final examination of the relations we have established among APE smearing, AUS 
smearing and cooling algorithms, we illustrate the topological charge density in Fig. 
Large positive (negative) winding densities are shaded red (blue) and symmetric isosurfaces 



TABLE III. The average of the ratio < napc(0.55)/napc(a) > or < naus(0.55)/naus(a) > for 
the action analysis, (5), and the topological charge analysis, {Q) from the 24^ x 36 lattice. 
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FIG. 20. Calibration of cooling and APE smearing via the ratio nc/raape(0-55) versus raape(0.55). 
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FIG. 21. Calibration of cooling and AUS smearing via the ratio nc/naus(0.55) versus naus(0.55). 
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aid in rendering the shapes of the densities. In fixing the x coordinate to a constant, a 
three-dimensional shce of a four- dimensional 24^ x 36 gauge field configuration is displayed. 
Fig. ^(a) illustrates the topological charge density after 21 APE smearing steps at a = 0.7. 
The other three quadrants display APE smearing, AUS smearing and cooling designed to 
reproduce Fig. |2^(a) according to the relations of Eqs. (p!7D, (|20|) and (pT]). The level of 
detail in the agreement is remarkable. 

V. CONCLUSION 

The APE smearing algorithm is now widely used in a variety of ways in lattice simula- 
tions. It is used to smear the spatial gauge-field links in studies of glueballs, hybrid mesons, 
the static quark potential, etc. It is used in constructing fat-link actions and in constructing 
improved operators with smooth transitions to the continuum limit. We have shown that 
to a good approximation the two-dimensionful parameter space of the number of smearing 
sweeps nape (a) and the smearing fraction a may be reduced to a single dimension via the 
constraint 

napc(Q^ ) _ , . 

^apc(a) a' ' 

satisfied for a and a' in the range 0.3 to 0.7. This result is in agreement with fat-link 
perturbation theory expectations, and survives for up to 200 sweeps over the lattice. We 
expect this relation to hold over the entire APE smearing range < a < 3/4. We find the 
same relation for AUS smearing provided a < 0.85. For AUS smearing, annealing of the 
links is included as one cycles through the Lorentz directions of the link variables in the 
smearing process. 

The smoothing of gauge field configurations is often a necessary step in extracting ob- 
servables in which the renormalization constants differ significantly from one, as is the case 
for the topological charge operator. It is also often used to gain insights into the nonper- 
turbative features of the field theory which give rise to the observed phenomena. We have 
determined that cooling, APE smearing and AUS smearing produce qualitatively similar 
smoothed gauge field configurations at /3 = 6.0 provided one calibrates the algorithms as 
follows: 

nc ^ 0.600 a ?2apc(«) , n,c ~ 0.618 anaus(tt) and anape(a) — 1.03a''riaus(a') • (26) 

The topological charge analysis serves to confirm the action analysis results at f3 = 6.0 
and further support these relations. However, it also reveals that at /? = 5.7 different trajec- 
tories are taken. At (3 = 5.7 different algorithms produce smoothed gauge field configurations 
with similar action, but different topological properties. 

As most modern simulations are performed at /3 > 6.0, the relations of Eqs. P3| ) and 
(^) will be most effective in reducing the exploration of the parameter space. It is now 
possible to arrive at optimal smearing by fixing the number of smearing sweeps and varying 
the smearing fraction, or vice versa. Finally, it is now possible to directly compare the 
physics of smeared and cooled gauge field configurations in a quantitative sense. 
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FIG. 22. Topological charge density of a 24^ x 36 lattice for fixed x coordinate. Positive 
(negative) windings are colored red to yellow (blue to green). Fig. (a) illustrates the topological 
charge density after 21 APE smearing steps at a = 0.7. Fig. (b) illustrates the topological charge 
density after 49 APE smearing steps at a = 0.3. Relation ( p!7[ ) suggests these configurations should 
be similar and we see they are to a remarkable level of detail. Fig. (c) illustrates the topological 
charge density after 20 AUS smearing steps at a = 0.7. Relation (p|) suggests these configurations 
should be similar and again the detail of the agreement is excellent. Finally (d) illustrates the 
topological charge density after 9 cooling sweeps, motivated by relation (]2^). While the level of 
agreement is not as precise, the qualitative features of the smoothed configurations are compatible. 
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We have used the action and topological charge trajectories to calibrate the various 
smoothing algorithms. While the topological charge may be related to the topological sus- 
ceptibility, it may be of future interest to introduce other physical quantities such as the 
string tension as a measure of smoothing ^j. 

It is well known that interpreting the physics of cooled configurations based on the 
standard Wilson action is somewhat of an art. The difficulty is that the 0(a^) errors in 
the action spoil an instanton under cooling by reducing the action below the one-instanton 
bound. A consequence of this is a lack of universality. For example, the action varies 
smoothly to zero while the topological sector jumps as (anti)instantons are destroyed. One 
may then attempt to define the concept of an optimal amount of smoothing, and perhaps 
also extrapolate back to zero smoothing steps. Unfortunately, all of these difficulties are 
apparent in the APE smearing analysis results as well. 

Because of the intimate connection between smearing and cooling, as emphasized in the 
discussion of Section |lllB 2 , it is natural to consider improved smearing, where the Hermitian 
conjugate of the improved staple is used in the smearing process. It will be interesting to 
see if the benefits of improved cooling carry over into improved smearing and this work is 
currently in progress. 
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